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INTRODUCTION 


A computer program based on state-of-the-art compressor ana 
structural technologies applied to bladed shrouded discs has been 
developed and made operational In NASTRAN Level 16. 

The problems encompassed Include aeroelastlc analyses, modes 
and flutter. 

The program is documented in the form of five NASA Contractor's 
Reports ■ ■ ■ one Technical Report and four Updates to NASTRAN Level 
16 Theoretical, User's, Programmer's and Demonstration manuals. 

This report describes the Theoretical manual updates. 


THEORETICAL MANUAL UPDATES 


18. AEROELASTIC , MODAL AND FLUTTER ANALYSES OF UNSTAllFO AXIAL FLOW TURBOMACNINES 


18.1 INTRODUCTION 

The rotors and stators of axial flow compressors and turbines are subjected 
to centrifugal, thermal and airloads that depend on the geometry and the operat- 
ing parameters. Steady aeroelastic and unsteady response of these "cyclically 
symmetric" structures, in turn, influence the applied thermal and airloads. 

These inter-active loads and responses arise fundamental ly from the elasticity 
of the structure and determine the performance and stability characteristics 
of the "flexible" turbomachin r 

Theoretical developments of References 1-3, have been applied to determine 
the thermal and airloads on the rotor/stator blade of an axial flow turbomachi ne . 
The computer code of Reference 1 has been adapted for NASTRAN in the functional 
module ALG to generate the steady state aerodynamic pressure and temperature 
loads. Computer codes of linearized, two-dimensional, harmonic cascade theories 
for subsonic and supersonic flows (References 2 and 3, respecti ve ly) have been 
utilized in the functional module AMG to estimate the harmonic airloads on the 
blade in a strip-theory manner. No transonic flow theory has been included 
presently, and the airloads on and near the transonic cylinder (or cone) are 
estimated by linear interpolation from subsonic and supersonic adjacent strip 
results. 

These steady and harmonic aerodynamic theories, in conjunction with the 
existing structural analyses capabilities in NASTRAN have been implemented in 
the form of two new rigid formats to perform; 

(1) Static aerothermoelastic "design/analysis", including differential 
stiffness effects, of an axial flow compressor rotor/stator 
(DISP Approach RF 16), and 

(?) Cyclic modal, unstalled flutter and subcritical roots analysts 

of an axial flow compressor and turbine rotor/stator (AER0 Approach 
Rf 9). 
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Th« rigid formats are designed such that the rotor (or stator) of a 
single-stage, or of each stage of a multi-stage compressor or turbine Is 
analyzed as an Isolated structure. 

Rigid formats have been designed In a modular fashion so that additional 
or improved aerodynamic computer codes could replace those currently 
incorporated. 
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18.? STATIC AEROTHERMOELASTIC ""ESIUM/ANAL YSIS" OF AXIAL 
FLOW COMPRESSORS WITH DIFFERENTIAL STIFFNESS 

At «n operating point under steady-stete conditions, the bladtd-disc 
of the compressor Is subjected to centrifugal, thermal an- aerodynamic 
loads that result In deformation of the elastic structure. For a fixed 
flow rate and rotational speed, the deformation implies a change In the 
operating point pressure ratio. 

The process of arriving at an "as manufactured" blade shape to produce 
a desired , design operating pcint pressure ratio at a given flow rate and 
speed is herein termed the "design" problem. The subsequent process of 
analyzing the performance ef the "as manufactured" geometry at off-design 
conditions ircluding the effects of flexibility is herein termed the 
"analysis" problem. 

The current NASTRAN Static Analysis with Differential Stiffness rigid 
format has been modified to include the effects of non-aerodynamic (centri- 
fugal, etc.) and aerodynamic (pressure and temperature) loads. The following 
remarks apply to the simplified problem flow and the algorithm shown in 
Fig urcs 1 and 2, respectively. 

1 The geometry of the compressor bladed-disc sector. Its Material pro- 
perties and the applied constraints are used to generate and partition 

the elastic stiffness matrix, Non-ae**odyanm1 c load vectors are formeo 
x . 

and an operating point flow rate, speed, loss parameters, etc. are selected, 

2. Based on the undeformed blade geometry and the operating point 
aerodyne mi c parameters, the functional module ALG generates the aerodynamic 
load vector. 

3. Total loads are defined as a combination of aerodynamic and non- 
aerodynamic loads. 

4. k linear solution for independent displacements is obtained based 
on the elastic stiffness and the total loads. 
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5. Omitted and constrained displacements are recovered, and stresses, 
reactions, etc., are obtained. 

6. A differential stiffness matrix is derived as a function of the 
gr id point displacements, 

7. A total stiffness matrix is now defined as a sum (or difference) 
of the elastic and geometric (differential) stiffness matrices for the 
"analysis" (or "design* 1 ) problem. 

8. The linear displacements obtained earlier are used to revise the 
blade geometry and a revised aerodynamic load vector is obtained. 

0. Again, the aerodynamic and non-aerodynamic load vectors are 
combined to define the total load vector. 

10. A non-linear solution for independent displacements is obtained based 
on the total stiffness and the total loads. 

11. Dependent displacements are obtained and data such as stresses, 
reactions, etc., are recovered. 

1 2. Convergence of the solution is based on the parameter c defined by 



Upon convergence, the final displacements, loads, he deformed blade 
geometry, etc., are output. Otherwise, further iterations are performed. 

A decision to update the differential stiffness matrix requires a shift 
to the outer loop. Only the load vectoi is revised in the inner loop iterations. 

12.1 The final pass, upon convergence, through the functional module 
ALC. yields tne flexible" operating point pressure ratio (among other 
aerodynamic data), which can be relocated on the compressor map. 
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The "design" node of the rigid format Is exercised only at the design 
operating point of the compressor. It is a two-step procedure in that 
having "designed" the blade shape, i.e., the "as Manufactured" shape. It 
should be "analyzed" at the same operating point to confirm the design point 
pressure ratio. The "analysis" mode of the rigid format is a one-step 
procedure. The "designed" blade is "analyzed” at selected operating points 
over the compressor map, one at a time, to generate the "flexible" performance 
characteristics of the compressor. The differential stiffness matrix 
generated during the analysis can be saved for use in subsequent modal 
analysis . 
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Differential stiffness matrix 1*^] 


Total stiffness matrix Uj*] 


Aerodynamic pressure and temperature 

loads (P *} on deformed blade, ALG 
9 


Total loads {P^J (Aerodynamic 
and non-aerodynamic) 


Ifdcpendent displacements {uj) 
(Non-linear Solution) 


Dependent displacements, stresses, etc, 
(Non-linear Solution) 


Adjustment to [K° g ]’ 
necessary 


.^Convergence 
Checks, DSCHK 


No change In £ 3 - 


Figure 1. Simplified Problem Flow for Static Aero thermoel as tl c 

"Des 1 gn/Ana lys i $* Rigid Format for Axial F lowCompressors 
including Differential Stiffness Effects, (continued) 
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1?.I 


1 ® 

Final displacements |u*), deformed 
blade geometry, strati, ate. ♦ 
operating point pressure ratio and 
othar flow parameters. 


cm 


I 


p 


FIgura I. Simplified Problem Flow for Static Aerotbermoblestlc 

■Design/Analysis* Algid Format for Axial Flow Compressors 
Including Differential Stiffness Effects, (concluded) 
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1 . 


2 . 


3. 


4 . 


Enter . after tha application of constraints and partitioning 
to tha stiffness Matrix and tha generation and transformation 
of the oea-aerodynamlc load vectors (centrifugal, etc.), with 



Aerodynamic Load Generator 
(pressure and temperature) 

<V * ♦(*/> 

{ p t ) ^ cojat r ai n { ^, 
partition 

* 

(U t ) • I *,,]" 1 


Undeformed blade geometry ♦ 
operating point (flow rate, 
speed, loss parameters, etc.) 


5 . (, a«» . »« . r . — (u t i 

[6^3 • [G 0 3. etc. 

6. It 4 "} +&!££* [kJ 9 ({u 9 ))3 
(Pg) • lPj A ) 


■0 


OUTER LOOP begins 


[K 


d 

aa 


] 


t constra In 
partition 



3 


u; ( 3 




{♦) for “analyst mode of the rigid format 
(-) for “design" mode of the rigid format 


<V ’ ‘V * i0> 

iuj) . 1u 9 ) 


Figure Simplified Solution Algorithm for Static Aerotheraot 1 as t i c 

“Design/Analysis* Rioid format for Arial F * ow Compressors including 
D i f f ert n t \ a i Stiffness Effects, (continued) 


•8.7-7 (9/30/78) 


AEROELASTIC AND DYNAMIC ANALYSES OF. TURBOMACHINES 


8. {P g ") 


■<s> 


-o 


B l lnner Loop begins 


Deformed blade geometry, revised with 
{u g *>. ♦ operating point. 


9. {P g2 > • {P g] > ♦ £P 9 -> 


{p b } ^constrain {p } 
1 partition 9 * 


10. {uj} -- [K^]' 1 (Pj > 


11. {u b }< ^recover 


^ etC - 


<•!> 


<«,*) ' C$1 


C, )■ C,) - c 5 ) 


t«$ g ) [«$, ((.$))) 


"ail* ’ c$) • rr, 0 » 


(l> „ 2 > ' <%„> * <V' 


12. Convergence checks ■« — ^ DSCH^ < {P g2 J» ^gI2^’ 


Differential Stiffness Checks 


E i V 


Exit with 


a. { u “ > , stresses, etc. 


Figure 2. Simplified Solution Algorithm for Static Aerothermoel asti c 

"Design/Analysis" Rigid Fermat for Axial Flow Compressors including 
Differential Stiffness Fffects. (continued) 
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b. 


Final deforced blade geometry 
+ operating point (flow rate, 
speed, loss parameters, etc.) 



{P 9 *) ♦ operating 


pressure ratio and 
other flow paraaeters. 


OR 2. 


e > e o 


and adjustment to K 


d 

99 


not necessary. 


Shift to the beginning of Inner Loop with 


*• <v ■ tf ,n* 


-0 


OR 3. c > c Q a_nd adjustment to K necessary. 


Shift to the beginning of Outer Loop with 


.. i. 5 i ■ t.,) 

b. 


-0 


Figure 2. Simplified Solution Algorithm for Static Aerotherir.oelastlc 

"Design/Analysis" Rigid Format for Axial Flow Compressors Including 
Differential Stiffness Effects, (concluded) 
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18.3 CYCLIC MODAL AND FLUTTER ANALYSES OF AXIAL 

FLOW TURBOMACHINES 

The problem of determining the complete, unstalled 
flutter boundaries of a cyclically symmetric compressor or 
turbine bladed disc involves each member set of the series of 
harmonic families of its modes, and the effects of permissible 
interblade phase angle, over an adequate set of operating 
points (flow rates, speeds, pressure ratios, implied Mach 
numbers, etc.) of the performance map. In view of the large 
number of variables influencing the definition of the flutter 
boundaries, a thorough parametric study requires systematic 
effective solution procedure. 

A capability, therefore, has been introduced in NASTRAN 
which, with repeated exercises over the range of variables 
involved, will enable determination of the flutter boundaries. 
The existing features of NASTRAN for Normal Modes Analysis using 
Cyclic Symmetry (Section 3,16, User's Manual) and Modal Flutter 
Analysis (Section 3.20, User's Manual) have been suitably 
combined for the cyclic modal, flutter and subcritical roots 
analyses in a new Rigid Format 9, Approach AER0. Provision is 
also made to include the differential stiffness effects by 
using the total stiffness matrix saved from the Static Aero- 
thermoelastic Analysis (see Section 18.2), 


In a compressor or turbine, an operating point 
Implies an equilibrium of flow properties such as density, 
velocity, Mach number, flow angle, etc., that vary across the 
blade span. Blade properties such as the blade angl e, stagger 
angle, chord, etc., also, in general, change from the blade 
root to the tip. The resulting spanwise variation in the local 
reduced frequency and the relative Mach number must be accounted 
for in estimating the chordwise generalized aerodynamic forces 
per unit span at each streamline. Integration of these forces 
over the blade span yields the blade generalized aerodynamic 
force matrix. Since the relative Mach number varies along the 
blade span, two two-dimensional, linearized, harmonic cascade 
theories (Refs. 2 and 3) one each for subsonic and supersonic 
flow have been implemented in a strip theory manner along the 
blade span. The chordwise aerodynamic matrices for streamlines 
with transonic inflow are derived by linear interpolation between 
those on adjacent (subsonic and supersonic) streamlines. 

The generation of the generalized air force matrices 
is an expensive operation and should be judiciously controlled. 

In the present development, the aerodynamic matrices are computed 
at a few reduced frequencies and interblade phase angles, and 
interpolated for others. Additionally, the chordwise generalized 

air force matrices are first computed for " aerodynamic modes " 
(heave, pitch, etc.). The matrices for chordwise structural 
modes are then determined from bilinear transformations along 
each streamline prior to the spanwise integration to obtain 
the complete blade generalized aerodynamic matrix. This permits 
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a change in the structural mode shapes of the same or a different 


harmonic number to be Included In the flutter analysis without 
having to recompute the modal aerodynamic matrices for aerodynamic 
modes . 


The following remarks apply to the simplified problem 
flow shown In Figure 1. In this figure, a compressor bladed 
disc performance map is shown, although the analysis Is equally 
applicable to both compressors and turbines. 

1. The geometry and the material properties of 

the bladed disc sector are defined along with the applicable 
constraints. An operating point is selected near the expected 
location of the flutter boundary. The solution procedure examines 
if this operating point is a flutter point. 

2. Flutter parameters such as densities, interblade 
phase angles and reduced frequencies are selected. 

3. The chosen operating point implies a certain 
spanwise variation of blade and flow properties. 

4. A harmonic number is selected for the cyclic modal 
analysis. Grid point mass and stiffness matrices are generated. 
The stiffness matrix saved from a previous Static Aerothermo- 
elastic Analysis can be used instead, and would include the 
differential stiffness effects at the steady state operating 
point under consideration. 


5. Constraints and partitioning yield the analysis 
set mass and stiffness matrices. 

6. Forward cyclic transformation results in the solu- 
tion set mass and stiffness matrices for the cyclic eigenvalue 
problem. 


7. Eigenvalues and eigenvectors in the solution set 
are obtained. 


8. Symmetric components eigenvectors are derived by 
a backward cyclic transformation. 

9. Symmetric components eigenvectors are augmented 
by recovering the dependent components, and are prepared for 
output if desired. 

10. For a non-zero harmonic number, the symmetric 
component eigenvectors are partitioned to separate the cosine 
and sine components. 

11. Based on the number of modes selected for flutter 
analysis the modal mass matrix is computed, 

12,13- Direct input mass, stiffness and damping matrices, 
if necessary, and the constraints thereon define these matrices 
for further analysis. 

14. The augmented eigenvectors, including any extra 
(or scalar) points introduced for dynamic analysis^are formed 
and used to define th^ new generalized mass, stiffness and 
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damping matrices. 


15. The streamline generalized aerodynamic matrices 
for chordwise aerodynamic modes are generated. The variation 
of the relative Mach number from streamline to streamline 
dictates the use of either of the subsonic and supersonic 
harmonic cascade theories. Such matrices for the streamlines 
with transonic inflow are interpolated. No transonic flow 
theory has been currently included. 

16. The structural modes are introduced via bilinear 
transformations along each streamline to define the chordwise 
generalized air force matrices. 

17. The blade generalized aerodynamic matrix is derived 
by a spanwise integration of the chordwise aerodynamic matrices 
for structural modes. 

18-20. The analysis loops through the user-selected 
combinations of density, interblade phase angle and reduced 
frequency. 

21. Based on the (o,k) combination, the appropriate 
blade aerodynamic matrix is chosen for the flutter equation. 
Linear or surface interpolation, at user's tion, is used if 
necessary . 

22. The generalized mass, stiffness and damping matrices 
of Step 14 and the generalized air force matrix of Step 21 are 
used to define the modal flutter equations. 



Y 

f 

y 

23. The solution to the flutter equations is 
sought In the form of complex eigenvalues and eigenvectors. 

24. The velocity-damping and velocity-frequency 
curves output for each (p,a,k) group are Interpreted to 
identify flutter points. 

25. Based on the relative stiffnesses of the blade and 
the hub of the bladed d 1 s sector, a series of harmonic numbers 
are investigated before arriving at the flutter boundaries. 
Presently, the solution rigid format is designed to accept 

one harmonic number at a time. 

The cyclic modal flutter analysis discussed herein 
is illustrated by the example 9-5-1 of the Demonstration 
Probl ems Manual . 
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